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The Hilbert-Chapman-Enskog expansion of the kinetic equations in mean flight 
times is believed to be asymptotic rather than convergent. It is therefore inadvis- 
able to use lower order results to simplify the current approximation as is done in 
the traditional Chapman-Enskog procedure, since that is an iterative method. By 
avoiding such recycling of lower order results, one obtains macroscopic equations 
that are asymptotically equivalent to the ones found in the Chapman-Enskog ap- 
proach. The new equations contain higher order terms that are discarded in the 
Chapman-Enskog method. These make a significant impact on the results for such 
problems as ultrasound propagation. In this paper, it is shown that these results 
turn out well with relatively little complication when the expansions are carried to 
second order in the mean free time, for the example of the relaxation or BGK model 
of kinetic theory. 
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I. INTRODUCTION 

Derivations of the dynamical equations for a continuum under the influence of internal 
friction were first produced in the nineteenth century by phenomenological arguments in 
the works of Navier, Poisson and Stokes. These "Navier-Stokes" equations contained co- 
efficients of viscosity and heat conduction whose evaluation was not provided by the early 
phenomenology. It was not until Maxwell's treatment of the problem through elementary 
kinetic theory that the dependences of these coefficients on the thermodynamical state vari- 
ables began to be understood. When the need for a kinetic theoretical underpinning of 
the theory of the dissipative properties of fluids had become clearer, Hilbert [1] developed 
approximate solutions of Boltzmann's basic kinetic equation in terms of an expansion in the 
mean flight time of the particles constituting the fluid. Finally, in the hands of Chapman 
and of Enskog, the derivation of the equations of fluid dynamics from the basic kinetic the- 
ory of gases took on the form that was most widely used throughout the last century. This 
theory, in the formal version that Enskog developed, was summarized authoritatively in the 
book of Chapman and Cowling [2], as well as by a number of later books [3-6]. 



*Electronic address: eas@astro.columbia.edu 
^Electronic address: jeanluc@mailaps.org 

^Present address: Department of Mathematics, Imperial College London, SW7 2AZ, United Kingdom. 



2 



Though it was satisfying that the same equations came from phenomenology as from 
kinetic theory, there was also bad news: when the mean flight time is long, the descriptions 
of sound propagation and of shock wave structure from the Navier-Stokes equations are 
not very good when compared to experiment [3]. This is one of the reasons that other 
approaches to the derivation of fluid dynamical theory have been sought over the years. The 
first approach was naturally to extend the Hilbert development to higher order. This was 
achieved by Burnett [7], who went to second and third order in the expansions. However 
the form of Burnett's equations that are conventionally employed are those that result when 
the Euler equations have been introduced into them as Chapman and Cowling [2] did (see 
the discussion by Agarwal et al. [8]). This may be the reason for the disappointing result 
that little or no improvement was found (and with often unphysical consequences) when 
the reduced Burnett equations were used for situations when the mean flight time was not 
extremely small. The reason that going to higher order does not lead to any real improvement 
in the results obtained with the Chapman-Enskog approach is probably that, as Grad [9] 
and Uhlenbeck and Ford [3] have suggested, the Hilbert expansion is asymptotic rather than 
convergent. Since the Chapman-Enskog procedure is a sort of iteration process, it would 
seem to be inappropriate for a series that is not convergent. That is, in the Chapman-Enskog 
procedure, one uses results from the previous order to simplify the results in the current 
order and this has the effect of restricting the domain of validity of the approximation at 
any order. 

To circumvent this limitation of the Chapman-Enskog procedure, an alternate approach 
was suggested in which the recycling of lower order results was avoided [10-13]. In the 
case of the relaxation model of kinetic theory [14, 15] this extension of the procedures was 
straightforward to apply and the results were encouraging. (The alternate approach was also 
carried out to first order for a Fokker-Planck type collision operator in Ref. [16].) Not only 
was an improved description of the structures of shock waves over that from the Navier- 
Stokes equations obtained, but also the phase speed of sound waves that emerged was in 
good agreement with experiment for all ratios of mean flight times to wave periods, unlike 
the Navier-Stokes results. Though the results for the damping of sound waves did not turn 
out as well, they were not bad for mean flight times below the sound period and they were 
better than the Navier-Stokes results. At present though, it is unclear whether a continuum 
theory can do much better than this result for damping. Perhaps a microscopic theory is 
required to do better, as Stubbe [17] and Sukhorukov and Stubbe [18, 19] suggest. Here we 
further probe the continuum theory by going to higher order using the new approach. 



We consider a gas of N identical particles of mass m whose individual positions and 
velocities are denoted by x and v. The distribution of the particles in the phase space with 
coordinates Xi and Vj = 1, 2, 3) is represented by a swarm of representative points whose 
local density /(x, v,t) is defined such that the probable number of particles in the volume 
element <ix dv is /(x, v, t) dx. dv. 

The fluid dynamical variables are defined in terms of moments of /. The density and 
fluid velocity are 



II. BASICS 




and 




(II.l) 
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and, with the peculiar velocity c := v — u, we may define the temperature 



3Rp 

We seek equations for these quantities. Interest has been focused on them because, as van 
Kampen [20] has stressed, these are the macroscopic "slow variables" that emerge from 
Boltzmann's theory. 

The particles obey Hamiltonian dynamics and so the phase fluid flows incompressibly 
and satisfies an equation of the form 

T>f = a-(3f (II.3) 

where 

V := d t + v% + a j d v3 (II.4) 

and summation over repeated symbols is understood throughout this work. The quantity a 
is the imposed external force per unit mass. 

On the left of (II. 3), we have the rate of change of / resulting from the streaming of the 
representative points through phase space. On the right, we see the effects of the interactions 
of the particles with those in their surroundings. Particles with velocity equal to v have 
their velocities changed to some other values at a rate f3f per unit volume while particles 
with other velocites are given the velocity v at the rate a. In strict equilibrium, / has the 
value /o, which does not depend on x or t. Then the left side of (II. 3) vanishes and we find 
that 

f = a/(3. (11.5) 

This expression recalls the Kirchhoff-Planck law of radiation theory and hints that radiative 
transfer theory is one source of inspiration for the relaxation model [21]. 

According to the description we adopt of the way the particles interact with each other, 
we may find / either as a consequence of the vanishing of the interaction term or we may 
need to prescribe it. In either case, the choice of greatest interest is the Maxwell-Boltzmann 
distribution 

n ( c 2 



/o(x,v,t) = exp 1^-—) , (II.6) 

where n := p/m. In the case where the fluid variables are allowed to depend on space 
and time, the use of the Maxwell-Boltzmann distribution for / defines a local thermo- 
dynamic equilbirium. When we introduce (II.5) (II.6) into the kinetic equation, with the 
understanding that / represents only a local equilibrium, we get the relaxation model 

Vf = ^°_zZ (II.7) 

where r := 1/(3. We may then study the case of nonequilibrium, with / and /o differing, 
by using the kinetic equation to follow the evolution of /. Alternatively, we may seek the 
coarser description of fluid dynamics in terms of the moments of the kinetic equation. 

Given the definitions of the fluid variables and the expression of fo in terms of them, 
we obtain consistency between the two by imposing the matching conditions [22, pp. 362 
and 369], 

r (fo-f)rd 3 v = (II.8) 



/ 
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where 

^ = m (1 , v l , \c 2 ) (II.9) 

with p — 0, 1, 2, 3, 4 and i = 1,2,3. This ensures that the state variables in / are the same 
ones that are given by the moments of /. Since we shall assume here that r does not depend 
explicitly on v, though it may depend on the state variables, the matching conditions are 
equivalent to the conservation of mass, momentum and energy. 

If we now multiply the kinetic equation by ^ and integrate over velocity space, we obtain 
the fluid equations 

d t p + V-(pu) =0, (11.10) 
<9 t u + u- Vu = -p -1 V-P + a, (11.11) 
\Rp{d t T + u- VT) = -P: Vu- V Q (11.12) 



where 

P 



:= Jmccfd 3 v and Q := J \mc 2 cf d 3 v (11.13) 

are the pressure tensor and heat flux vector. The trace of the pressure tensor is 3p where 

p:=RpT; (11.14) 

here, R = k/m is the gas constant and k is Boltzmann's constant. 
Another useful quantity is 

S:=C v \n(P-) (11.15) 



where C v is the specific heat at constant volume and 7 is the ratio of specific heats; 7 = 5/3 
in this work and C v = The expression for S is that for the specific entropy of a perfect 
gas when p is the pressure. On using the continuity equation and the equation of state, we 
see that 

S 



(d t + u-V)[ — ) =AlnT+fV-u. (11.16) 



We may then rewrite the heat equation (11.12) as 

pT (d t S + u • VS) = - (P - pi) : Vu - V • Q (11.17) 

where I is the unit tensor. This result motivates the introduction of a notation for the tensors 
of kinetic theory, namely that an x over any tensor designates its trace while a circle over 
a tensor, say A, signifies its symmetrized and traceless part, as in 

A := |(A + A) - |AI, (11.18) 
where the overtilde indicates a transpose. For example, we see from the definition of the 

X o 

stress tensor that P = 3p and P = 0. 

To obtain determinate expressions for P and Q, we seek an approximation for /. In what 
follows we describe a procedure for obtaining and using one. We begin with a standard 
expansion of / in terms of r: 

/ ./o • /1 r ■ f, t 2 ■ .... (11.19) 
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When we introduce f n T n (no summation) into the formulae for P and Q we obtain contri- 
butions written as P( n ) and Q*") which lead to 

P = p(°) + p( 1 ) + ... and Q = Q (0) + Q {1) + .... (11.20) 

At the lowest order, in which f is the Maxwell-Boltzmann distribution, we find that 
p(°) = p i an d Q(o) = o, and so we have (11.10) and 

<9 t u + u- Vu + p~ 1 Vp-a= 0(r), (11.21) 
lRp(d t T + u- VT)+pV-u= 0(r). (11.22) 

Without the (r) terms, these are the Euler equations. 

At the next order, we obtain the fluid equations found earlier [10, 11]; these reduce to 
the standard Navier-Stokes equations when r becomes strictly infinitesimal. As these first- 
order developments are central to the current work we present them next (Section III). This 
provides us with an opportunity to describe how the Navier-Stokes equations emerge and to 
indicate how the Lorentz force enters the theory when electromagnetic effects are considered. 
Then we proceed to carry the development on to the next order (Section IV) and to show the 
predictions of the second-order theory for the propagation of free sound waves (Section V). 

III. FIRST-ORDER THEORY 

A. Approximation for / 

At order r, we readily obtain from Eq. (II. 7) that 

A = -£/„■ (III.l) 
It proves advantageous to multiply (III.l) by l// and to rewrite it as 

A = -Dln/ , (IH.2) 
where /i := /i//o- Then, with (II. 6), we have 

A = -s> 



c 2 , n 
+ In 



2RT (2%RTfl 2 



(III.3) 



which we may write out as 



h^c-Vc-C^-^T-^,. (III.4) 



B. Developments 

Some further rearrangements are usefully made on noticing that 

Dc = D(v-u) = a-Au-(c-V)u (III.5) 
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where D t — d t + u • V. We can then rewrite Eq. (III. 4) as 



/i = -^c-[a-Au-(c-V)u 

c 2 



- (AhiT + c- VlnT) -D t \np-c- Vlnp. (III.6) 



2RT 2 / 

It is useful at this point to introduce the rate-of-strain tensor 

E:=|(Vu + Vu) (III.7) 

X 

where the tilde over a tensor indicates transposition. Then, the trace E = V ■ u and we may 
rewrite the c • (c • V)u term on the right of (III. 6) as 

o i 

c-(c- V)u = cc : E+ -c 2 V • u. (III.8) 

After introducing the definition of S and using the ideal gas law (11.14) and the continuity 
equation (11.10), we find by a longish series of straightforward substitutions that (III. 6) 
becomes 

f l = -( VlnT--i-cc:E + c-*- ( - ^] £ (III.9) 
J \2RT 2) RT \2RT 2 J y J 

where 

E:=A(|) (HI.10) 

and 

*:=-^ (Au-a)-Vlnp. (III.ll) 



On using (11.11), we may also write 



* = iv-P (111.12) 

V 



In a final cosmetic touch, we define the dimensionless peculiar velocity 



w := c/V2RT (111.13) 

and rewrite (III. 9) as 

/i = - (w 2 - - j c- VlnT-2i«;^ : E + c • * - ( w 2 - - J E . (111.14) 



C. The Prandtl Number 



The relaxation model is prone to some ambiguity coming from the the lack of specificity 
in the details of the interactions of particles. In particular, the model leaves open the way r 
depends on v. This is really a matter that is related to the nature of the atomic interactions 
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and does not seem to have a significant qualitative influence on the macroscopic results. 
However, the Prandtl number, a, given by the relaxation model is different from that pre- 
dicted by the Boltzmann model which is restricted to two-body interactions. (The Prandtl 
number is the ratio of the viscosity times C p to the thermal conductivity.) Experiments in 
noble gases give values for a close to the value 2/3 predicted by the Boltzmann equation for 
a gas of hard spheres; by contrast, the value unity is found with the relaxation model. In 
calculating macroscopic quantities with the relaxation model, notably the stress tensor and 
energy flux, it is sometimes considered desirable to adjust the value of a in the results to be 
2/3. Since there is more than one way to do this, a systematic procedure seems called for, 
particularly when the Lorentz force comes into play. 

Woods [23] has given a prescription for introducing the desired value of the Prandtl 
number that seems to work well, though it is more of a rule of thumb than a fully justified 
procedure. It has the advantage of providing a definite prescription that can be carried into 
higher order and it leads to agreement with results from the Boltzmann equation for all of 
the numerical coefficients that appear in the higher moments obtained from the Chapman- 
Enskog procedure in second order in r. Since this modification of the results can be retracted 
by setting a — 1, it is harmless and we shall introduce it. 

In Woods' approach one introduces the Prandtl number into the formulae for approxi- 
mating / rather than into the final fluid equations. The motivation for Woods' procedure 
is that the times of relaxation for momentum and energy are different [24] so that he has 
suggested the use of different relaxation times depending on whether a term from / enters 
the pressure tensor or the heat flux. Terms with even powers in c contribute to the pressure 
tensor and are assigned a relaxation time r; terms with odd powers in c contribute to the 
heat flux and are assigned a relaxation time r/a where a is the Prandtl number. Therefore, 
we rewrite (III. 14) as 

A = ~ (w 2 - ^ c- VlnT- 2ww : E+ c • (w 2 - S. (111.15) 

This way of adjusting the Prandtl number differs from that in the first-order derivation 
of Refs. [10, 11] and, for a = 2/3, it leads to agreement of our second-order theory with 
Burnett's results [7] from the Boltzmann equation if we discard terms of order r 3 (see the 
appendix). 

If the acceleration a is not a function of v, Eq. (III. 15) is cubic in d, so f\ is a polynomial 
in c\ If a is a polynomial in v, then we can still write f\ as a polynomial in c, possibly of 
higher order than cubic. An interesting instance of this arises when a includes the Lorentz 
force on charged particles. The only term in (III. 15) that involves a is That term vanishes 
for the standard Chapman-Enskog expansion without the electric field E and magnetic field 
H. When E and H are included, ^ is nonzero and can depend on velocity through a 
[Eq. (III.ll)]. Let 

a(x, v, t) = — [E(x, t) + v x H(x, *)] (111.16) 

m 

be the acceleration due to the Lorentz force, where e is the electron charge in MKS units. 
By rewriting a as 

a(x, v, t) = — [E(x, t) + u x H(x, ()] + -cx H(x, t) . (111.17) 

m m 

we split the acceleration into a part that is independent of c and a part linear in c. As we 
see from Eq. (III. 15), a will be dotted with c, and so c x H will not contribute. Hence, 
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electric and magnetic fields can be included by simply expressing the acceleration as 



which is independent of v. 



a(x, t) = — [E(x, t) + u x H(x, t)} , (111.18) 

m 



D. Matching Conditions 



Although the relaxation model is often regarded as linear, in much the same way that the 
analogous equation of radiative transfer is, it has to be recognized that fo depends on the 
moments of the full / and so cannot be considered as given a priori. Treating the problem 
as formally linear is possible because the feedback on the fluid variables from alterations in 
/ are forced to be consistent with the basic equations by the matching conditions (II. 8). In 
the present work, they guarantee the conservation of the collisional invariants ip^. 

If we insert the solution (III. 15) for j\ into (II. 8) we find that the density matching 
condition is satisfied to this order. The momentum matching condition is that 

mv /id 3 ?; = ^^* = - 1 -V-(P-pI). (111.19) 

2(7 2(7 

But P(°) = pi so that P — pi = (r). Hence, the momentum matching condition is satisfied 
to the order of the current approximation, which is all that can be asked in an asymptotic 
development. 

Similarly, the energy matching condition is that 

\ J mv 2 /i d 3 v = -\pRT E. (111.20) 

with E given by (111.10). From (III.ll) and (11.12), we find that (111.20) is 

E = D t InT + | V • u = -— P : Vu + V • Q . (111.21) 

6 3p 

On recalling that = 0, we see that E is (r). We thus conclude that all the matching 
conditions are satisfied to the required order. 



E. Pressure Tensor and Heat Flux 

As mentioned in Section III, if the acceleration a does not depend on v (except possibly 
for a magnetic term) and r does not depend on v, then we can write /i as a third degree 
polynomial in the components of v . Finding the contributions of the first-order terms in / 
to the pressure tensor and heat flux is then just a matter of inserting the solution (III. 15) 
for fi (times r) into (11.13) and carrying out the resulting integrals, as indicated by (11.20). 
We find 

pW = -pr(2E + El), Q (1) = -f ^-pRT {VlnT - V) . (111.22) 
Here we see the appearance of the viscosity and the conductivity, 



jj,:=pr and rj:——C p pT 

<7 



(111.23) 
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where we have written C p for 5R/2. The Prandtl number is a = fiC p /i]. 

The first terms in each of the expressions (III. 22) are those that appear in Navier-Stokes 
equations. The second terms (proportional to S and *S?) are the terms that are retained when 
the additional approximations introduced in the usual Chapman-Enskog procedure are not 
made. That is, if we had used the zeroth-order results (the Euler equations) to eliminate the 
time derivatives contained in £ and we would have obtained the Navier-Stokes equations, 
as in the Chapman-Enskog method, though generally speaking this extra (and unnecessary) 
step is normally taken earlier in the development in the standard approach. Our suggestion 
is that if one wants to simplify the expressions for P* 1 ) and QM by eliminating the time 
derivatives, it is more accurate to use the current approximation for the equations of motion 
rather than a lower-order one. If this elimination is done iteratively, it shows that our 
expressions for the higher moments contain terms of all orders in r. However, we have left 
the time derivatives as they appear in £ and \l/ [Eqs. (III. 10) and (III. 11)], and we shall 
simply work with the full system self-consistently when we illustrate a physical problem in 
Section V. 

The stress tensor and energy flux at this stage are written as 



where P (0) = pi and Q<°) = 0. Since the terms by which and differ from the 
corresponding Navier-Stokes expressions are 0(r 2 ), the differences are within the allowed 
error estimates for both treatments. By the same token, we find directly from the definition 
of P that its trace is 3p while the trace of our expression for it differs from 3p by terms of 
order r 2 . This is also compatible with the asymptotics to this order. As we shall see, in the 
next order, the trace of our approximation to P differs from 3p only by terms of order r 3 . 
A suggestive way look at these results is to observe that we may write 



To this approximation, if we replace p by the first two terms in its Taylor series in time, we 
have 



That is, the form is like that of the Navier-Stokes equations except that the pressure is 
evaluated at one relaxation time prior to the present, and there is a volume (or bulk) 
viscosity term. 



p = p(0) + p(l) + ( r 2) Q = Q(0) + Q(l) + o ( r 2) 



(111.24) 




(111.25) 




(111.26) 



IV. SECOND-ORDER THEORY 



A. 



The Distribution Function 




(IV.l) 



which in terms of / 2 := /2//0 becomes 



/ 2 = -£>/!- /if In /o -ADmr. 



(IV.2) 
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We allow for a dependence of r on the macroscopic state variables, p and T, but assume 
that r does not depend on the particle velocities, either microscopic or macroscopic. We 
then write 

Vlnr = rD t \np + sD t \nT (IV.3) 



where 



<91nr 



and 



<91nr 



(IV.4) 



dlnp dlnT 

Chapman and Cowling specify r = — 1 on the authority of Maxwell while, depending on the 
density of the gas, s varies from about 0.5 to 1 [2, p. 223]. 

With r = — 1, the third term on the right-hand side of (IV. 2) is 



/iSlnr = /i (sDlnT-Dlnp). 
After introduction of the fluid equations (11.10)— (11.12), this becomes 

/i V lnr = A [|(5 - 2s) V • u + (s - 1)2 + s c ■ V InT - c • V hip] 
The middle term on the right of Eq. (IV. 2) may be written out as 



(IV.5) 



(IV.6) 



/ 1 Dln/ = - 



w 2 - I I c • V InT + 2w w : E - - c • * + I w' z - 



x 



w 2 c- VlnT + 2ww :E-c-*+ [ w 2 - 



(IV.7) 



The first of the three terms on the right of Eq. (IV. 2) requires a certain attentiveness 
so it is well to divide its calculation into manageable pieces. Firstly, the derivative of the 
peculiar velocity (III. 5) can be written as 



Dc = a 2 (* + V hip) — c • E + c x f2, 



where 



a 



2 ._ 



2RT, 



and O := |V x u is the vorticity. When we observe that 

Va = \ot (E — | V • u + c • V In T) ; 

we find that (IV. 8) implies 

= V (-^—) = a(* + Vlnp) -w-E + w x ft - \ w{Z + c • V InT). 

Because it is antisymmetric, the vorticity term disappears from Dc 2 : 
T)c 2 = 2c-Vc = a 2 (c ■ * + c V lnp) - 2c c : E , 

and similarly 

T>w 2 = a(w- * + w ■ Vlnp) -2ww : E - w 2 (S + c • V InT) . 



(IV.8) 
(IV.9) 

(IV. 10) 

(iv.il) 

(IV.12) 
(IV. 13) 
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Next we have 

T 2) V In T = D t VT + c • VVT + § V • u VT - V In T (c • V)T — EVT. (IV. 14) 

The derivatives of E and are then expanded in the usual manner for v-independent 
quantities as 

2)S = AS + c-VE, = A* + c ■ V*. (IV.15) 

The first term on the right-hand side of (IV. 2) becomes 

1 
D/i = (a(w • * + u> • Vlnp) - 2™w : E - w 2 (£ + c • VlnT))c • VlnT 

a 

~a ( U ' 2_ ^) (l ft2 (* + Vln ^) ~ C ' E + C x n ) • VlnT 

-~ ( w2_ C ' (AVT + c- VVT + |V-uVT- VlnT(c- V)T-SVT) 

- 4w ^±a(* + V lnp) -w-E + wxQ-± u>(£ + c • V lnT)^ : E 

- 2u> u; : (£> t E + c • VE) + - (|a 2 (* + Vlnp)-c-E + cx n) ■* 
+ - c • (A* + c • V*) - (w 2 - I) (AS + c-VE) 



- (aw • (* + Vlnp) — 2tu it> : E — w 2 (£ + c • V lnT))£ . (IV.16) 

Solving for f 2 is now a matter of collecting the terms of (IV. 7), (IV.6) and (IV.16). We 
then follow the same prescription as in Section III and assign a relaxation time r to the terms 
in f 2 that are even in c, and r/er to the odd terms, to account for the different relaxation 
rates of momentum and energy. Once this is done, we may verify, in an analogous manner to 
the procedure of Section III D, that the matching conditions are satisfied by f 2 to the order 
of the approximation. This is a lengthy, but straightforward computation that we omit here. 



B. Pressure Tensor and Heat Flux 



As in Section HIE, the second-order pressure tensor and heat flux are obtained by in- 
serting the solution f 2 of Section IV A into (11.13). On performing the required integrations 
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we find, for the pressure tensor, 



O 

( ( ° \ a 2 s a 2 

F {2) =pr 2 \ ^A(2E + EI)-2Vu-(2E + EI)J +— VVT + — VTVT 

o 

° ° 4 ° 5o;^ 5^ o? 1 4 ° ° 

+ 8E ■ E + -(f - s) V ■ uE + — -V 2 TI + — — | VT| 2 1 + - E : EI 
3 00" 1 ooi ^ 3 



9 o 9 o 9 o 9 

a „ s . — or ° — a or 



+ 2(s-5)EE + — (17 -3s)* VlnT+ — * Vlnp+ — ** V* 

2o o o o 



a 2 „ a 2 



+ sE 2 I+ — * Vlnpl+ — (| VlnTI 



3o o 
1 

3 



1 rv 2 W "4 

-(5 -2s) V-uSI+^^ 2 I-9-V-*I , (IV. 17) 
3 3o oo J 



and for the heat flux, 

Q( 2 ) = ipaVo" 1 (5(T- 1 -(AVT-Vu-VT) + 14 (s + (1 + a' 1 )) E ■ V InT 

VT ° ° 

+ f a ~ 1( y\ - s)V -u^ + 14 V - E-4E - Vlnp 

- 5o" _1 (A* -Vu *) + 10VE - 2(2 + 7o" _1 ) * ■ E 
- 5(1 + (s + l)o- x ) E * - o^f (5 - 2s) V ■ u E * 

+ 5((f + 2s) + so" 1 ) EVlnT — 5EVlnpj . (IV.18) 

A circle on top of a bar means that the entire term under the bar should be symmetrized 
and rendered traceless as implied by the notation of (11.18). 

We may expand the second-order pressure-tensor (IV. 17) and heat flux (IV.18) in r, 
to recover the same terms as one obtains using the Chapman-Enskog approach, that is, 
what is usually called the BGK-Burnett equations. In doing this, we discard terms of 
order r 3 , include the second-order contributions from and and treat the Prandtl 
number in the same way in both cases. When we put in a Prandtl number of 2/3, and 
discard cubic terms, we recover the conventional Burnett result for the Boltzmann collision 
operator, which is a justification for the Prandtl number fix discussed in Section III C. To 
illustrate the nature of the new form for the pressure tensor and heat flux, we treat the case 
of ultrasonic sound waves in the next section. 



V. ULTRASONIC WAVES 

It is known experimentally [25] that, when the frequency of a sound wave becomes compa- 
rable with the collision frequency of particles, both the Navier-Stokes and Burnett equations 
give poor accounts of the dispersion relations. In Refs. [10, 13, 26] the present theory was 
found to give very good agreement with the data for phase velocities, in first order, so that 
it is not possible to do much better with higher-order terms to within experimental error. 
However, the first-order theory predicts no damping of the waves at large Knudsen num- 
ber. One of our motivations for this work was to see if we could capture the experimental 
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damping in higher order. As we shall see, the second-order theory does better than the first 
order theory for a Knudsen number up to about 8. The price for this seems to be negative 
damping for Knudsen numbers above 8. This instability of the homogenous state appears 
also in the Burnett theory and it remains to be seen whether any of the regularizations 
proposed for that case will work for our equations. 

We start by expanding the macroscopic variables about an equilibrium rest state (u = 0), 

p = p o (l + 0), T = T (l + 6), p = p Q (l + m), (V.l) 

and find from the equation of state (11.14) that w = <fi + 9. The linearized equations of 
motion (II.10)-(II.12) are then 

d t (f) + V • u = 0, (V.2) 
Po d t u + V • (P (0) + P (1) + P (2) ) = 0, (V.3) 
l Po d t 9 + Po V • u + V • (Q« + Q( 2 )) = 0, (V.4) 

where the linearized pressure tensor in each order is 

F {0) =p wl; P {1) = -/i (Vu+ Vu + d t 9l); (V.5) 



P {2) = Po?o [(l + a- 1 )d t (Vu + Vu) + -9 t V-uI 
v a 

2 2 

+ Q 2 9 1 + ^ ( VV0 + I V 2 1) + ^ (VVzu + \ V 2 w I) ] , (V.6) 

and the linearized heat flux is 

Q(D = - Vo T V(0 + zu)-lfi a' 1 d t u ; (V.7) 

Q( 2 ) = 1^ a 2 r a' 1 (-^d 2 u + 5(2 + a- 1 )d t V9 + -d t Vw + 7V 2 u + 9V(V-u)) . (V.8) 
v aa a ) 

The viscosity fi and the thermal conductivity rj are defined by Eq. (III. 23). The derivatives 
of t do not enter the linearization, as we see from the absence of the parameter s defined 
in (IV.4). 

We take the divergence of the u equation (V.3), and use the continuity equation (V.2) to 
eliminate V • u. After some manipulation we find 



rl $4 = a 2 r 2 V 2 (0 + 9) + \ (2«9 t V 2 - d t V 2 9) 

CL 

4 ( -(2 + 3a' 1 ) d 2 V 2 (f> + 8 2 V 2 9 + Qa 2 ^ 1 V A 9 + SaV" 1 V 4 ^] , (V.9) 



r 4 

+ 



where 



2 PO PO f\r -\n\ 

a := — , u := — (V10) 
Po Po 

are respectively the isothermal sound speed and kinematic viscosity, in terms of 
which r = Po/Po — Vo/a 2 . The equation for the temperature perturbation (V.4) becomes 

| r d t 9 - t d t <t> - |a 2 r 2 a' 1 V 2 (29 + </>) + § r 2 a' 1 d 2 <P 

+ \Tla 2 a- 1 [ %d^ + 10(l + a- 1 )d t V 2 9 + (5a- 1 -16)d t V 2 (j)) =0. (V.ll) 
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If we use t as unit of time and ar as unit of length, but do not change the meaning of 
the symbols used otherwise, we can rewrite (V.9) and (V.ll) as 

(<9 2 - V 2 - 2d t V 2 + (2 + 3a- 1 ) <9 2 V 2 - 3a- 1 V 4 ) 0+(<9 t -l-<9 2 -6 ( x- 1 V 2 )V 2 fl = 0, (V.12) 

and 

(3d t - lOa-V 2 + 10a- 2 (a + 1) 8 t V 2 ) 9 

+ 0~— 1 (5<9 2 - 5V 2 - 2ad t - Sa^d? + <r _1 (5 - 16a) d t V 2 ) = 0. (V.13) 

Taking the density and temperature perturbations and 9 proportional to exp(i(kx — out)) 
we obtain a relation between k and uj, assuming neither of and 9 vanishes identically. The 
number e = ujtq is dimensionless and can be regarded as the Knudsen number, the ratio of 
the relaxation time (tq) to the macroscopic time The approximation involved in a 

Chapman-Enskog-type expansion get worse as ut becomes large. 

To compare to experimental results, we plot the inverse phase speed k/uo as a function 
of 1/utq. Figure 1 shows the real part, and Fig. 2 the imaginary part (damping). The 
phase speed is in units of the adiabatic sound speed, (tpo/ Po) 1 ^ 2 , where 7 = 5/3. We 
used a Prandtl number of 2/3 and so compare with the N-S and Burnett results from the 
Boltzmann equation for Maxwell molecules. 

VI. DISCUSSION 

According to Grad [9] and Uhlenbeck [3], the Hilbert expansion of / in terms of mean 
flight times is asymptotic and not convergent. For this reason, we have avoided using an 
iterative procedure like those of Chapman and of Enskog. In their ways of deriving the fluid 
equations [2], it is the practice to use lower-order approximations to simplify the current 
approximation. In particular, the Chapman-Enskog approach uses the Euler equations to 
simplify the calculation of the fluid equations at first order, and the Navier-Stokes equations 
at second order. Here, we do not recycle lower-order results in this way and so we obtain 
equations at a given order of accuracy iV that differ from the Chapman-Enskog results 
by terms (r Ar+1 ). In other words, our results agree asymptotically with those of the 
Chapman-Enskog procedure but, when we apply the two theories at the edges of their 
validities, we see differences in their quantitative predictions. 

The effect of this difference in approach is seen starkly when we take the trace of the 
stress tensor, P. From the definition of the pressure tensor, its trace is seen to be 

P = 3p (VI. 1) 

while our iVth approximation to the stress tensor, F N , has a trace of 3p + 0(t n+1 ). It 
appears to be a common misconception that it is required to make the 0(t n+1 ) error in the 
trace of the pressure tensor vanish identically. However, we are perfectly at liberty to retain 
terms 0(r N+1 ) in a theory that is intended to be good to 0(r N ). Indeed, we suggest that, 
in the present instance, it might even be desirable to retain the extra terms in P. 

In the approach described here, we have in the order N, an approximation for the pressure 
tensor with an error of order r N+1 , so that we may write P = Pat + 0(r N+1 ). Let us put 
this as 

F = pI + Ar N + 0{r N+1 ) (VI.2) 
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where A is the tensor we need to find to complete our approximation. When we take the 
trace of this equation, we find the requirement that 

A + O(r) = 0. (VI.3) 

Iterative procedures like Chapman-Enskog theory force the trace of A to vanish exactly 
and, on that account, (VI.3) cannot be satisfied unless r is truly infinitesimal. But, in our 
approach, (VI.3) is not an embarrassment since it points a way to increase the domain of 
validity of our approximation. By making a suitable choice of A we might cancel some of the 
remaining error in the trace of P if we but knew how to arrange this. As we have observed, 
the expansion procedure offers a natural choice for A and the problem is then to test its 
accuracy. 

As in most physics problems, we prefer to let the choice among the theories be made by the 
experiments, which do favor our approach over the standard Chapman-Enskog treatment. 
Indeed, with the Prandtl number adjustment, our results from the relaxation model seem to 
be preferred to those from the Boltzmann equation in the Chapman-Enskog approach. We 
shall return to the discussion of our approach to the Boltzmann equation in another place. 
Here, we have remained in the relatively simple case of the relaxation model whose results 
seem as good as those with the Boltzmann model once the Prandtl number is adjusted (see 
also the discussion of the so-called BGK-Burnett equations in [8]). 

In the study of ultrasound, the propagation speed we predict agrees well with the experi- 
mental findings at all Knudsen numbers. However, our results for the damping of ultrasound 
are good (as judged by experiment) only for Knudsen number up to order unity. Since our 
second-order results are an improvement over the first order ones, we may hope for im- 
provement as higher orders are achieved. But in practice, going to higher order than second 
would seem to be prohibitively complicated, even if the calculations could be performed. 
The problem is then what ought to be the approach to the damping of sound waves in the 
limit of very long mean flight times of particles. 

One clue is that for moderate Knudsen numbers we have some negative damping for a 
small range of flight times. This is reminiscent of the instability found in the Burnett equa- 
tions and Jin and Slemrod [27, 28] treat this problem with a regularization procedure that 
involves manipulation of high order terms, as Rosenau has also done [29] . Such procedures 
have also proved effective in problems of photon transport [30]. Similar procedures may be 
used to ameliorate our equations and we shall return to that possibility elsewhere. But there 
is an interesting aspect to this. The various regularization procedures being introduced to 
improve the continuum equations generally involve the introduction of higher order terms. 
Here, we are making a plea for first looking more carefully at the higher order terms that 
are discarded by standard procedures. However, we wish to conclude with mention of a final 
issue. 

Stubbe [17] has argued that the dynamics of rarefied gases cannot be well treated macro- 
scopically and a that microscopic approach is needed as in the case of plasma physics where 
no macroscopic approach has produced Landau damping [19, 31]. On the other hand, it 
may be that a suitable resummation of the results of lower orders would resolve the issue. In 
any case, it appears that the study of the damping of ultrasound waves is going to involve 
serious complications. 
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FIG. 1: The real part of the inverse phase speed, k/uo, versus 1/lotq. The phase speed is given in 
units of the adiabatic sound speed, (tpo/po) 1 ^ 2 - The first- and second-order dispersion relations of 
this paper (with a = 2/3) are compared to Navier-Stokes, Burnett, and the experiment of Meyer 
and Sessler [25]. 
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APPENDIX A: RECOVERING BURNETT 

As we have said, the Chapman- Enskog procedure eliminates terms of higher order than 
does the current procedure. This elimination can be performed at various stages and, when 
we eliminate those same terms from our equations, we recover the Chapman-Enskog type of 
results in the N-S and Burnett levels of approximation. More surprsingly perhaps is that 
when we set the Prandtl number to 2/3 the same elimination of higher order terms gives 
us the Burnett equations with exactly the same coefficients as is found with the Boltzmann 
equation for hard spheres. Here, we sketch this reduction. 

Equations (IV.17)-(IV.18) contain terms in E and which are formally of order r. We 
may thus write 

o o 
( ( o ox /t,2 o q,2 ° o o 

P (2) = pr 2 2 I D t E — 2Vu • E I + — VVT + — - VTVT + 8E-E 
^ v ) aT aT 2 

+ -(|- S )V.uE+^V 2 TI+^|VT| 2 I + -E:Ell + (r 3 ) (A.l) 

and 

Q( 2 ) = ipaVff" 1 (5<7 _1 -(AVr-Vu-VT) + 14 (s + (1 + a" 1 )) E • V InT 

+ f ( 7- 1 (|-s)V-n-^ + 14V-Ej +0(r 3 ). (A.2) 

These are not the second-order pressure and heat flux of the Burnett equation. To recover 
the Burnett result, we need to expand to first order in r the £ and \l/ terms appearing in 
the first order pressure tensor and heat flux (III. 22). 

First we do this for St, which enters the heat flux in (III. 12). On inserting (III. 22) 
in (III. 12) we find 

* = I v (-2p r E) + (t 2 ) (A.3) 
p 

which can be expanded to give 

^ = _2r V • E - 2rE- Vlogp- 2rpE • Vlogr + (r 2 ) . (A.4) 

These terms, when multiplied by (5r/4o")po; 2 and added to (A.2), agree with the Burnett 
expression for the second-order heat flux derived using a Boltzmann collision operator (after 
setting a = 2/3). 

Now for E, which enters the pressure tensor in (III. 12). After inserting (III. 22) in (III. 21) 
we obtain 

E = - — (-2prE)V -u + V • (-\-pa 2 V\nT) + (r 2 ) ; (A.5) 
3p V 4 a / 

Upon multiplying this by —prl, adding to (A.l), and setting a = 2/3, we recover the 
Burnett expression for the second-order pressure tensor appropriate for the Boltzmann col- 
lision operator. Thus our full expressions for the pressure tensor and heat flux agree with 
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Burnett to second-order, as they must since both theories (ours and Burnett's) are second- 
order expansions. The perfect agreement also validates the Prandtl number fix discussed in 
Section III C. 
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